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ABSTRACT 

Aims. We propose the Wind of Fast Rotating Massive Stars scenario to explain the origin of the abundance anomalies observed in 
globular clusters. 

Methods. We compute and present models of fast rotating stars with initial masses between 20 and 120 M Q for an initial metallicity 
Z=0.0005 ([Fe/H] a -1.5). We discuss the nucleosynthesis in the H-burning core of these objects and present the chemical composi- 
tion of their ejecta. We consider the impact of uncertainties in the relevant nuclear reaction rates. 

Results. Fast rotating stars reach the critical velocity at the beginning of their evolution and remain near the critical limit during the 
rest of the main sequence and part of the He-burning phase. As a consequence they lose large amounts of material through a mechan- 
ical wind which probably leads to the formation of a slow outflowing disk. The material in this slow wind is enriched in H-burning 
products and presents abundance patterns similar to the chemical anomalies observed in globular cluster stars. In particular, the C, N, 
O, Na and Li variations are well reproduced by our model. However the rate of the 24 Mg(p, y) has to be increased by a factor 1000 
around 50 x 10 6 K in order to reproduce the whole amplitude of the observed Mg-Al anticorrelation. We discuss how the long-lived 
low-mass stars currently observed in globular clusters could have formed out of the slow wind material ejected by massive stars. 

Key words. Nuclear reactions, nucleosynthesis, abundances - Stars: rotation - Stars: mass-loss - Stars: abundances - globular clusters: 
general - globular clusters: individual: NGC 6752 



1. Introduction 

Galactic globular clusters (hereafter GCs) appear to be chemi- 
cally homogeneous (with the notable exception of co Cen) with 
respect to the iron-group (Mn, Fe, Ni, Cu), neut ron-capture 
(Ba, La, Eu) and alpha -el ements (Si, Ca) (e.g., Kraft et all 
Il992t Uames et alJl2004albt ISnedenll200l ISobeck et al.ll2006l) . 
However it has long been known that these large aggregates 
of stars show strong inhomogeneities in lighter elements: C, 
N, O, Na, Mg and Al abundances present indeed large star-to- 
star abundance variations within all the individual GCs studied 
up to now (for complete references see t he early rev iews by 
iFreeman & Norrisl 1 19811 ISmithlll987l a nd iKraft] 1 1 994l and th e 
more recent ones bv lGratton et al.l 20041 and ICharbonnell2005l) . 

The observed patterns point to the simultaneous operation of 
the CNO, NeNa and MgAl cycles of hydrogen-burning : C and 
N, O and Na, and Mg and Al are respectively anticorrelated, the 
abundances of C, O and Mg being depleted while those of N, 
Na and Al are enhanced. Whenever C, N, and O are observed 
simultaneously, their su m appears to be con s tant within the ob - 
servational errors (e.g., iDickens et all 119911: llvans et al.|[l999l) . 
The sum Mg+Al is also found to be constant in several clus- 
ters dShetrond [19961) . Observations in NGC 6752, M 13 and 
M 7 1 show that the Mg depletion is due to the burning of 24 Mg 
while 25 Mg is untouched and 26 Mg is produced in the Al-rich 
stars (lYong et al.1120031 120051 l2006h . Last but not least, Li was 



found to be anticorre lated with Na in turn-off stars of NGC 6752 
(IPasquini et al.ll20"05t) . All these features are considered anoma- 
lous because they are not seen in field stars of similar metallicity 
(e.g jGratton et al.ll2000l) . 

For many years only the brightest GC red giants were ac- 
cessible for detailed spectroscopic observations, and two main 
theoretical streams were competing to explain the available data 
: (1) the so-called "evolution" scenario according to which the 
chemical anomalies are generated inside the low-mass stars we 
are presently observing, and (2) the "self-enrichment" (or pri- 
mordial) scenario according to which such patterns pre-existed 
in the protocluster gas and were inherited at the birth of the long- 
lived stars. 

The evolution hypothesis has been seriously challenged by 
recent spectroscopic observations of less luminous stars in ear- 
lier stages of evolution in a number of GCs. Such studies re- 
vealed indeed that stars located slightly above and below the 
main sequence tur noff exhibit the same anomalies as their 
giant counterparts (iGratton et al. 20011: Grundahl et al. 20021: 
Carretta et al.l(2003l 120041: ICohen et al.ll2002tlRarmrez & Cohen! 



20021 2003; Ha rbeck et a!]|2003l) . However such objects are not 



hot enoughQ for the required set of nuclear reactions to occur 
within their interior. Let us recall that while the CNO cycle is ac- 
tivated for temperatures above 20 x 10 6 K, the NeNa and MgAl 
chains require temperatures around 35 x 10 6 K and 50 x 10 6 K 
respectively. Destruction of 24 Mg by proton-capture needs still 
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1 In the central region of a 0.85 M e , [Fe/H] = -1.3 turnoff star, the 
temperature is of the order of 25 x 10 6 K. 
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highe r temperatures, around 70 x 10 6 K (e.g., lArnould et al.l 
1 19991: Prantzos & Charbonnel in preparation). As a consequence 
the abundance variations cannot be produced in situ, but cer- 
tainly reflect the initial composition of the protostars. It is thus 
clear now that a large fraction of GC low-mass stars were formed 
from material processed through H-burning at high tempera- 
tures and then lost by more massive and faster evolving stars, 
and perhaps mixed with some original gas. Various aspects 
of this "self-enrichment scen ario" are discussed in details by 
iPrantzos & Charbonnell d2006l hereafter PC06). 

Regarding the nucleosynthetic site, most studies have fo- 
cused on massive AGB stars which were sugg ested as the pos- 
sible polluters by ICottrell & Da Costal (11981 1) . The two main 
reasons why these objects have been favored are that (1) they 
host regions where H-burning occurs at high temperatures (in 
particular when they experience the so-called hot bottom burn- 
ing, or HBB, at the base of their convective envelope be- 
tween successive thermal pulses), and (2) the material they 
eventually eject (by stellar winds or Roche lobe overflow) 
is not enriched in iron. This last property is well consistent 
with the observational fact recalled above that the iron abun- 
dance in a GC does not show any significant scatter. For long 
the AGB hypothesis was discussed only on a qualitative ba- 
sis. However recent custom-made stellar models (Ventura et al. 



2001, 2002; Denissenkov & Herwig 2003; Karakas & Lattanzio 
20031: lHerwig|l2004allbl: IVentura & D'Antondl^OOlallbllcl l2006t 
Decressin et al., in preparation) pointed out very severe diffi- 
culties from the nucleosynthesis point of view which stem from 
the competition between the HBB and the third dredge-up. This 
latest process does indeed contaminate the envelope of the AGB 
with the products of helium burning and crea tes abundance pat- 
tern s in conflict with th e observed ones (see iFenner et al.ll2004l 
and ICharbonnelll2005l for more details). PC06 discuss other 
shortcomings of the AGB scenario, the main one being related 
to the peculiar initial mass function it requires. In addition, they 
underline the fact that the AGB scenario gives no satisfactory an- 
swer as to the role of stars more massive and less massive than 
the presumed polluters. 

Massive stars, m ore precisely, Wolf- R ayet s tars, have 
been proposed by Brown & Wallersteinl (1 19931) and by 
IWallerstein et al.l (1 19871) as possible sources for the very 
early enrichment of g lobular clusters. More recently 
iMaeder & Mevnea (l2006f) suggested that He-rich stars in 
a> Cen could be formed from wind materi al of fast rotating 
massive stars. IPrantzos & Charbonnell (120061) proposed a com- 
prehensive (albeit qualitative) scenario for the role of massive 
stars, suggesting that their winds provide the metal-enriched 
material for the next stellar generation, and that the subsequent 
supernova explosions provide the trigger for the star formation; 
the SN ejecta escape altogether the GC environment, along 
the cavities opened previously by the stellat winds. PC06 also 
studied the massive star IMF required to explain quantitatively 
the amount of Na-enhanced stars observed in NGC 2708 and 
found it to be flatter than canonical (i.e. Salpeter) IMFs; even 
flatter IMFs would be required in case the polluters wer e AGB 
stars. Similar conclusions for a flat IMF are reached by ISmithl 
(120061) . for the case of N enhancement of GC by massive star 
winds. The main reason why such objects have been discarded 
in the past in the context of the self-enrichment scenario is 
related to the fact that iron is ejected at the time of their super- 
nova explosion and this constitutes a priori a serious drawback 
for considering massive stars as GC pollution sources. This is 
true unless some filtering process removes the ejecta enriched 
by helium burning and more advanced nuclear stages while 



preserving those bearing the signatures of hydrogen processing. 
This is the key point of the Winds of Fast Rotating Massive Stars 
scenario (hereafter WFRMS) that we propose in the present 
work. In our framework, the GC chemical anomalies are built in 
H-burning zones of massive stars. Rotational mixing brings to 
the surface CNO-processed material which can then be ejected 
in a slow wind when the stars rotate a t the cr itical limit. 

As sho wn bvlSackmann & Anandl(ll970l) and later bv lLangerl 
d 1998b and lMaeder & Meynej d200lb . massive stars do reach the 
so-called critical velocitjo early on the main sequence if (1) they 
start their evolution with a sufficiently high initial rotation rate, 
(2) they do not lose too much angular momentum through stellar 
winds and (3) an efficient mechanism (meridional circulation in 
our models) transports angular momentum from the core to the 
envelope. Once the critical limit is reached, the surface veloc- 
ity remains near the critical value during the rest of the main 
sequence and very likely an equatorial disk is formed as ob- 
served for instance around Be stars. In the case of these stars, 
only minor o utflow in the line-formin g region is observed (see 
the review bv lPorter & Riviniusll2003l) . This material has thus a 
great chance to be retained in the GC potential well. On the other 
hand fast rotation leads to strong internal mixing of the chemi- 
cals. As a result the ejected material will present the marks of H- 
processing occurring in the stellar core. Thus if new stars form 
out of the slow wind material, their composition would bear the 
signatures of H-burning. The main question addressed in this 
paper is whether the chemical composition of the ejecta of fast 
rotating massive stars is compatible with what observed in the 
long-lived GC stars. 

In the present paper we develop in details the WFRMS sce- 
nario. In § 2 we describe the physical ingredients of our models 
of rotating massive stars. We first focus on the properties of the 
60 M Q models computed with various assumptions. For this ini- 
tial stellar mass we discuss the nucleosynthesis in the H-burning 
core in § 3 while in § 4 we investigate the chemical composi- 
tion of the wind ejecta. The full range of masses between 20 
and 120 M is investigated in § 5. In § 6 comparisons between 
the wind composition and the GC abundance patterns are per- 
formed. A schematic and speculative discussion of the complete 
scenario for explaining the inhomogeneities in GCs is presented 
in § 7. The conclusions and some future lines of research are 
given in § 8. 



2. Physical inputs 

Our stellar models were computed with the Geneva evolutio n 
code including the effects of rotation dMevnet & Maederll2000b . 
We focus on the mass range corresponding to stars with high 
enough central temperatures on the main sequence for the NeNa 
and MgAl chains (see Fig.[TJ to be activated. Models with initial 
masses of 20, 40, 60, and 120 M are computed from the Zero 
Age Main Sequence up to the end of the core He-burning phase, 
and a 200 M Q model is computed up to the end of the core H- 
burning phase. The theoretical predictions for the 60 M star are 
discussed in detail in § 3 and 4 in order to properly illustrate the 
WFRMS scenario. 



2.1. Microphysics 

We use the OPAL opacities dlglesias & Rogers! I1996I) . 
complemented at temperatures below 5000 K with the 



2 By critical velocity, we mean the equatorial surface velocity such 
that the centrifugal acceleration exactly balances the gravity. 
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molecular opacities of lAlexander & Fergusonl 
(http : //webs Twichita . ed u/physics/opacrEyj). 




d 19941) Table 2. Nuclear reaction rates adopted for the NeNa- and 
MgAl-chains in the s ets B-D. \N] [H 2] and [H041 corre- 
spo nd respectively to lAngulo et al] dl999l) . llliadis et aiT(l200ll) 
and lHale et ail d2002l (20041) . "Norn.", "low." and "up." refer re- 
spectively to nominal, lower and upper limits of the experimen- 
tal values. For all the other reactions of our network we use the 
NACRE nominal values 



Fig. 1. Network of nuclear reactions involved in the NeNa and 
MgAl chains. Unstable nuclei are in shaded boxes. Arrows indi- 
cate the type of nuclear reactions: (p, y), f3 + , and (p, a) 



Table 1. Main physical inputs of the stellar models for the var- 
ious initial masses considered. The labels "r" and "s" indicate 
respectively the models computed with or without rotation while 
the labels A-D refer to the choices in nuclear reactions (see text 
and Table 0). The initial value of Q/ Q cr j t is given . [I], [ H02] 
and [H041 correspond respectively to llliadis et all d200ll) and 
lHaleetalid2002l 120041) 



Reaction 


set A 


set B 


set C 


set D 


i(, Ne(p, y) 


[N], nom. 


[N], nom. 


[N], low. 


[N], low. 


21 Ne(p, y) 


[N], nom. 


[I], nom. 


[I], low. 


[I], low. 


22 Ne(p, 7 ) 


[N], nom. 


[H2], nom. 


[H2], low. 


[H2], low. 


23 Na(p, 7 ) 


[N], nom. 


[H4], nom. 


[H4], low. 


[H4], low. 


23 Na(p, a) 


[N], nom. 


[H4], nom. 


[H4], up. 


[H4], up. 


24 Mg(p,7) 


[N], nom. 


[I], nom. 


[I], up. 


[I], +3 dex 


25 Mg(p,7) 


[N], nom. 


[I], nom. 


[I], up. 


[I], up. 


26 Mg(p, T ) 


[N], nom. 


[I], nom. 


[I], up. 


[I], up. 


27 Al(p, 7 ) 


[N], nom. 


[I], nom. 


[I], low. 


[I], low. 


27 Al(p, or) 


[N], nom. 


[I], nom. 


[I], low. 


[I], low. 



while in set C some specific rates are set to the experimental up- 
per or lower limits. Fig. [2] presents the corresponding rates for 
the temperature range between 30 x 10 6 and 80 x 10 6 K which 
is typical of the central temperatures of our main sequence stars. 
Finally set D is similar to set C except for the proton-capture on 
24 Mg which is i ncrea sed by three orders of magnitude compared 
to llliadis et aD(l200lh nominal value around 50 x 10 6 K. 

The initial composition of the chemical mixture is given in 
Table [3] It corresponds to tha t used to compute the opacity ta- 
bles (Igl esias & Rogers! [19961 Weiss alpha-enhanced elements 
mixture). The metallicity of our models is [Fe/H] ^ -1.5 corre- 
sponding to that of NGC 6752 which is the GC with the largest 
set of abundance data. The initial isotopic ratios of magnesium 
are taken equal to 80:10:10; this corresponds to the values ob- 
served in NGC 6752 "unpolluted" stars (i.e., in stars with high O 
and low Na abundances) in contrast with "polluted" stars which 
display larg e O depletion with high Na abundance dYong et al.l 
120031 120061) . 



M(M ) 


Label 




Nuclear rates 


60 


60rA 


0.95 


set A T 


NACRE (nominal) 




60rB 


0.95 


setB': 


[I, H02, H04] (nominal) 




60rC 


0.95 


setC + : 


[I, H02, H04] (exp. limits) 




60rD 


0.95 


set D f 






60rE 


0.80 


setC 




20 


20rC 


0.95 


setC 




40 


40rC 


0.98 


set C 




120 


120rC 


0.80 


set C 




200 


200rC 


0.95 


set C 




20 


20sC 





set C 




40 


40sC 





set C 




60 


40sC 





set C 




120 


120sC 





set C 





Table 3. Initial abundances in mass fraction 



Element Abundance Element Abundance 



l H 


0.754 




1.53e-8 


3 He 


2.93e-5 


20 Ne 


5.30e-5 


4 He 


0.245 


21 Ne 


5.00e-8 


12 C 


3.50e-5 


22 Ne 


4.72e-6 


'- ! c 


1.47e-7 


23 Na 


3.30e-7 


I4 N 


1.03e-5 


24 Mg 


1.68e-5 




1.58e-8 


25 Mg 


2.10e-6 


16Q 


3.00e-4 


26 Mg 


2.10e-6 


l7 o 


1.31e-7 


27 A1 


9.00e-7 


18Q 


7.45e-6 


28 Si 


2.56e-5 



Details on the nuclear rates used are presented in Table[2] 



In order to investigate the effects of the nuclear reaction rate 
uncertainties on the products of hydrogen nucleosynthesis, we 
present four models for the 60 M star computed using different 
sets of nuclear reactions for the hydrogen-burning network (see 
Table[T|). Set A uses all the n ominal values of the NACRE com- 
pilation (Anguloet al. 1999). The t hree other cases include the 
updates of llliadis etatl d2001l) and lHale et alj d2002l |2004[) for 
the reactions involved in the NeNa and MgAl chains but with 
different options (see Table |2|. Set B includes nominal values 



2.2. Rotation and mass loss 

We fo llow the formalism by IZahnl d 19921) and iMaeder & Zahnl 
dl998l) for the transport of angular momentum and chemicals 
in rotating stars. The effects of both meridional circulation and 
shear turbulence are taken into account: the meridional circu- 
lation advects angular momentum and the shear acts as a dif- 
fusive process. The transport of chemical species is computed 
as a diffusive process as the result of meridional circulation and 
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Fig. 2. Nuclear reaction rates for NeNa (left) and MgAl chains (right) from the set A, B and C (top to bottom). Arrows and upper 
full line in the lower right panel indicate the increase of the rate of 24 Mg(p,y) 25 Al at T — 50 and 60 x 10 6 K assumed in set D. 



horizontal and vertical turbulence dChabover&Zahnll9 92). The 
treatment of the convective instability is done according to the 
Schwarzschild criterion and we do not consider overshooting. 

The tr eatment of rotation inclu des the hydrostatic effects 
following iMevnet & MaedeJ d 19971) as wel l as the impact of 
rotatio n on the mass loss rate described by iMaeder & Mevneil 
d2000t) . We d o not account fo r the wind anisotropics induced by 
rotation as in IMaeder! (1 19991) although the related effects would 
reinforce the trends found in this paper by fastening the arrival 
at the break-up limit (see § 4. 1). 

The radiative mass loss rates are from Kudritzki & Pulsl 
(120001) when log T eff > 3.95 and from lde Jager et al.l (1 1988b oth- 
erwise. When a model reaches the WR phase (i.e., when the sur- 
face hydrogen mass fraction becomes lower than 0.4 and the ef- 
fective temperature is higher t han 10 4 K), the mass lo ss rate is 
switched to the prescription of iNugis & Lamersl d2000l) . Except 
for the WR phase, we consider a dependence of the mass loss 
rates with metallicity asMoc VZ/Z , where Z is the mass frac- 
tion of h eavy elements at the surface of the star. 

As in lMevnet et alj d2006l) a specific treatment for mass loss 
has been applied a t brea k-up. Let us recall that according to 
IMevnet & Maederi d2000t) three kinds of "break-up limits" can 
be defined depending on which mechanism (i.e., radiative ac- 
celeration, centrifugal acceleration or both) contributes to coun- 
terbalance the gravity: 1 .- The r-Limit, when radiation effects 
largely dominate; 2- The Q-Limit, when rotation effects are es- 
sentially determining break-up; 3.- The QT-Limit, when both 
rotation and radiation are important for the critical velocity. In 
the present work the O-Limit is reached during the MS, while 
the QF-Limit is encountered by our most massive stellar models 
(M > 60 Mo) just after the MS. 

During the MS, near the Q-Limit, two counteracting effects 
are competing. On one hand, matter is removed from the stel- 
lar surface (mainly through the equatorial regions) together with 
angular momentum. Also the expansion of the envelope tends to 
slow down the surface. On the other hand, meridional advection 



in the outer layers acts so as to transfer angular momentum from 
th e inner stellar regions to the surface (see for instance Fig. 1 
in IMevnet & Maeder 2002). This acts as to accelerate the sur- 
face. As long as the internal transport of angular momentum is 
efficient enough for accelerating the outer layers in a timescale 
shorter than the mass loss or inflation timescale, the surface ve- 
locity remains near the critical limit. In practice, however, the 
critical limit contains mathematical singularities; we thus con- 
sider that during the break-up phase, the mass loss rate is such 
that the rotation velocity stays near a constant fraction of the 
critical value (0.98 typically). At the end of the main sequence, 
the stellar radius inflates so rapidly that meridional circulation 
cannot anymore ensure the internal coupling and the break-up 
phase ceases naturally. However the most massive stellar mod- 
els encounter the QF-Limit after the MS. In that case, we ap- 
ply similar procedures as the one described above, i.e. we apply 
mass loss rates that maintain the model at a constant distance 
from the Or-Limit. When the star has lost sufficient amount of 
mass for evolving away from this limit, the regular mass loss 
rates apply again. 

For the reasons invoked in the introduction and which will 
be discussed again later, we explore the case of high initial rota- 
tional velocities. For the stars with initial mass between 40 and 
120 M Q we take V ini = 800 km s _1 . The 200 M Q model starts 
with y^i = 1000 km s _1 . This corresponds to an initial value 
of Q/Ocrit between 0.80 and 0.98. The 20 M star starts with 
Vini = 600 km s _1 as it nearly corresponds to the break-up ve- 
locity (Q/fioi, = 0.95). 

3. Central hydrogen-burning in a 60 M star 

As a prerequisite to our scenario we first have to check whether 
the abundance patterns due to nuclear reactions in the hydrogen- 
burning core of massive stars do mimic the chemical trends ob- 
served in GC low-mass stars. Obviously the abundance varia- 
tions obtained within the central regions are the most extreme 
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Fig. 3. Evolution of the central abundances in the 60 M rotating 
models computed with different sets of nuclear reaction rates: 
set A (full line), set B (long dashed lines), set C (short dashed 
lines) and set D (dotted lines). The amount of the unstable nu- 
cleus 26 Al is added to those of 26 Mg and Mg. 



ones that one can expect within the WFRMS scenario. Indeed 
some dilution of this processed material is expected in the ra- 
diative envelope of the polluters (see § 4) and then eventually 
later with the intracluster gas (see § 6). We are thus looking 
for stronger abundance variations in the stellar core than the 
presently observed ones for the elements involved in the CNO, 
NeNa and MgAl cycles. 

We first investigate the nucleosynthesis that occurs within 
the convective core of a main sequence 60 M Q star and we ex- 
plore in detail the uncertainties of the nuclear reaction rates. The 
evolution of the central abundances of the key elements is pre- 



sented in Fig. [3] as a function of the central abundance of hydro- 
gen for the various sets of nuclear reactions (see Table|2]i. 



3. 1 . CNO nucleosynthesis 

In the 60 M Q star the central temperature on the main sequence 
varies from 48 X 10 6 K to 75 x 10 6 K (see Fig.|7). The CNO cycle 
thus rapidly reaches equilibrium at the beginning of H-burning. 
As a result the central abundance of 12 C drops suddenly by about 
an order of magnitude, while that of 14 N increases by a factor of 
29. Slightly later 16 does reach its equilibrium value, which is 
1 .6 dex below its initial abundance. After this adjustment phase 
at the very beginning of the main sequence, the CNO elements 
stay at their equilibrium level. These predictions are identical for 
all the 60 M models presented here because they were all com- 
puted with NACRE prescriptions for the CNO-cycle reaction 
rates (no modification of these rates has been published since 
the NACRE compilation was made available). 



3.2. NeNa and MgAl predictions with the NACRE nominal 
values - Model 60rA 

Model 60rA was computed with the NACRE nominal values for 
the whole nuclear network (full lines on Fig. [3]). 

As can be seen on Fig. [3] the central abundance of 23 Na 
shows a three steps evolution on the main sequence: a first rapid 
raise by 1 .2 dex, then a more progressive increase of 0.7 dex until 
the mass fraction of hydrogen at the center, X c , equals 0.4, and 
finally a decrease. The initial feature is due to proton-captures on 
21 Ne and 22 Ne which are the fastest reactions of the NeNa chain 
(see Fig. |2j». As proton-capture on 23 Na is slower, the abundance 
of sodium increases until the complete consumption of 21 Ne and 
22 Ne. Later on 20 Ne starts burning in favor of sodium. The slow 
increase is due to the competition between the burning of 20 Ne 
and that of 23 Na via the channels (p,y) or (p,a). When using 
NACRE nominal values the transformation of 20 Ne into 23 Na 
is more efficient than the 23 Na destruction for T inferior to ~ 
50 x 10 6 K which results in a slow increase of 23 Na. When the 
temperature exceeds 50 x 10 6 K the situation reverses and the 
abundance of 23 Na decreases. 

Let us now turn our attention to Mg and Al. In model 60rA 
the central abundance of 24 Mg increases slightly during the main 
sequence. That of 25 Mg first decreases by nearly 4 dex, then it 
presents a sawtooth behavior until X c equals ~ 0.2 and finally 
it increases. 26 Mg first slightly increases before decreasing dur- 
ing most of the core H-burning phase and increasing again when 
X c becomes lower than ~ 0.07. The behavior of Mg (total) fol- 
lows mainly that of 24 Mg. Last but not least, the abundance of Al 
increases. 

The increase of the 24 Mg abundance is due to proton- 
captures on 23 Na. The very efficient destruction of 25 Mg is due 
to the reaction 25 Mg(p,y) 26 Al. The sawtooth behavior results 
from small instabilities affecting the size of the convective coreQ 
The slight production of this element at the end of the main se- 
quence is due to the proton-captures on 24 Mg. The initial raise of 
the abundance of 26 Mg results from the destruction of 25 Mg de- 



3 When the size of the convective core slightly increases, some 25 Mg 
is dredged down. For most elements, the quantity dredged down is small 
compared to the abundance in the core; however in the case of 25 Mg 
which is severely depleted in the core, this temporarily produces a small 
increase of the central abundance. 
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scribed aboveQ When all the 25 Mg is consumed within the core, 
26 Mg is converted into aluminum. 

Table |4] sums-up the mean variations of the central abun- 
dances during the main sequence for the NACRE nominal re- 
action rates. O and Na are respectively depleted and produced 
as required by the data. On the other hand the increase of the 
total magnesium we obtain is at odds with the o bserved Mg-Al 
antico rrelation. These results agree with those of lArnould et al.l 
(119991) who investigated in details hydrogen-nucleosynthesis at 
constant temperature and solar metallicity with NACRE reac- 
tion rates. They found a production of Na and they argue that it 
is dominated by uncertainties on the rates for T > 50 x 10 6 K. 
Also 24 Mg is not destroyed even when experimental errors on 
rates are taken into account for T < 70 x 10 6 K. Production of 
aluminum is due to the burning of both 25 Mg and 26 Mg. 

3.3. NeNa and MgAI predictions with the llliadis nominal 
values - Model 60rB 

In model 60rB we take into account the nomin al rates given by 
llliadis et al l d200ll) and lHale et all J20021I2004 for the reactions 
involved in the NeNa and MgAI chains. The differences with the 
NACRE prescriptions can be seen in Fig. [2] At T = 50 x 10 6 K 
the rates of the 22 Ne(p, y) and 23 Na(p, y) reactions are lowered 
respectively by 1.7 dex and 1.3 dex compared to NACRE nom- 
inal values, while the rate of 23 Na(p, a) is 0.6 dex higher. For 
27 Al(p, y) the new rate is increased by 1 dex at T — 50 x 10 6 K, 
but is unchanged for T > 70 x 10 6 K. Finally the rate of the 
27 Al(p, a) 24 Mg reaction is now higher by 0.7 to 1.1 dex in the 
range of temperature of interest. The prescriptions for the other 
reactions are unchanged. 

This update does not affect the global structure and evolution 
of the star, as the energizing rates concern the CNO cycle which 
is not modified. Thus the lifetime, the size of the convective core 
as well as the mass loss history are identical in all our models. 

The key reaction for the NeNa chain is 20 Ne(/?, y) as 20 Ne is 
the most abundant neon isotope and its destruction by proton- 
captures is much slower than that of the 2I Ne and 22 Ne (by 2.7 
and 3.7 orders of magnitude respectively). Since the rate of this 
reaction is unchanged the behavior of the Ne isotopes is barely 
affected with respect to the results presented in § 3.2. On the 
other hand because the rate of proton-capture on 23 Na is reduced 
now with respect to the former case, the consumption of sodium 
is lowered (its abundance does not decrease at the end of the 
main sequence). As a consequence the abundances of 24 Mg and 
of the total magnesium do not increase although one does not ob- 
tain the decrease required by the observations. Finally the 26 Mg 
burning is favored so that aluminum is produced earlier on the 
main sequence although its final abundance reaches the same 
value as in model 60rA. 

In summary the new set of nuclear reactions resolves partly 
the difficulties encountered by model 60rA regarding the anticor- 
relation between Mg and Al. However the predicted magnesium 
isotop i c ratio s are still in conflict with the data by lYong et all 
(l2003ll2005l) . 

3.4. NeNa and MgAI predictions with experimental extreme 
values - Model 60rC 

In order to try and solve the remaining difficulties, we look now 
for the most favorable set of nuclear reactions by considering the 

4 The abundance of the unstable isotope 26 Al which decays into 26 Mg 
is included in the total 26 Mg abundance. 



published experimental limits. In model 60rC (set C) we make 
the following assumptions (see Table |2}: We take the lower lim- 
its for the burning rates of all the neon isotopes in order to lower 
the overall production of sodium. We also take the lower limit 
for the 23 Na(p, y) reaction which becomes now the slowest re- 
action of the NeNa chain. This reduces the linkage between the 
two chains and disfavors the production of 24 Mg. Regarding the 
MgAI chains, we take the upper limits for the destruction rates 
of the magnesium isotopes and the lower limits for the burning 
rates of 27 Al. 

The corresponding predictions are shown by the short- 
dashed lines in Fig. [3] The overall increase of the 23 Na abun- 
dance is more modest than in the previous models. One ob- 
tains now an important decrease of the 26 Mg abundance, while 
the final predictions for 25 Mg and 27 Al are essentially not af- 
fected with respect to sets A-B. Again, 24 Mg stays constant dur- 
ing the hydrogen-burning phase. Nevertheless the total magne- 
sium abundance decreases by 0.1 dex in this case. The stronger 
variation observed could come from more massive stars where 
the central temperature is hotter and the burning rate of 24 Mg is 
faster (see § 5). 

3.5. Model 60rD 

As shown previously the difficulty regarding the Mg destruc- 
tion comes mainly from the fact that the central temperature 
reaches the required extreme values only at the very end of the 
main sequence (see Fig. |7). We thus note that the use of the 
rate published for the 24 Mg(p, y) reaction iPowell et al.l d 19991) 
is not compatible with significant Mg depletion within the core 
of massive main sequence stars. This led us to tentatively modify 
this rate in order to reconcile the theoretical predictions with the 
abundance data. 

In model 60rD we keep all the reaction rates as in 60rC ex- 
cept for the burning rate of 24 Mg that we artificially enhance by 
a factor of 10 3 and 10 LS at 50 and 60 x 10 6 K respectively (the 
recommended rate is used for temperature lower than 40 x 10 6 K 
and higher than 70 x 10 6 K). With these assumptions this reac- 
tion becomes as efficient as the other ones involved in the MgAI 
chain (see the arrows on Fig.|2]i. 

As a result (see the dotted lines in Fig. |2]i the 24 Mg abun- 
dance decreases by about 0.4 dex when the central hydrogen 
mass fraction equals ~ 0.2. At the end of the main sequence 
25 Mg and 26 Mg are respectively slightly destroyed and produced. 
Simultaneously the total magnesium abundance decreases by 
0.3 dex while that of 27 Al increases by 1 .2 dex for a central hy- 
drogen mass of 0.2. These predictions are in good agreement 
with the observational constraints. 

3.6. Summary 

Our models were computed with the metallicity of NGC 6752 
([Fe/H]=-1.5) with the aim of comparing our theoretical predic- 
tions with the observations in this GC. In this section we have 
focused on the nucleosynthesis within the core of a 60 M main 
sequence star. The corresponding abundance variations are thus 
the most extreme ones one might expect, and they cannot be 
compared directly with the observational data. However we see 
already that we hold a serious candidate for building up the ob- 
served inhomogeneities and that it is worthwhile to explore fur- 
ther the WFRMS scenario. 

The situation can be summarized by a close inspection of 
Table|4]where we list the most extreme abundance variations be- 
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Table 4. Mean core abundance variations during MS of the 
60 M models. 5X refers to the amplitude (in dex) of variation 
of element X. For NGC 6752 it refers to the amplitude variation 
observed between polluted and unpolluted stars 



Model 


SC 


<5N 


60 


<SNa 


SMg 


6A\ 


NGC 6752 


-0.7 


1.7 


-1.0 


0.9 


-0.3 


1.4 


60rA 


-0.9 


1.5 


-1.8 


1.6 


0.5 


0.7 


60rB 


-0.9 


1.5 


-1.8 


1.6 


0.0 


0.7 


60rC 


-0.9 


1.5 


-1.8 


1.0 


-0.1 


0.7 


60rD 


-0.9 


1.5 


-1.8 


1.0 


-0.4 


1.3 



tween polluted and unpolluted stars in NGC 6752 as well as in 
the core of our 60 M models. We see that the comparison is 
rather favorable as far as the O-Na anticorrelation is concerned 
and that sets A and B lead to higher Na production than set C. 
On the other hand, the Mg-Al anticorrelation appears to be more 
difficult to reproduce. The use of sets A and B lead to a Mg-Al 
correlation, while that of set C builds up a weak Mg-Al anticor- 
relation. The situation becomes however more favorable when 
one increases the rate of the 24 Mg(p, y) 25 Al reaction as in model 
60rD. In that case 24 Mg is destroyed while the abundance of 
25 Mg is barely affected and that of 26 Mg slightly increases. This 
is in agreement with the observationa l constraints obtai ned on 
the Mg isotopes in NGC 6752 stars bv lYong et ai1d2005l) . 

Let us note that the H-burning signatures we describe are 
similar to those found in the H-burning shell during central He- 
burning. Indeed in this shell the temperature does not surpass 
60 x 10 6 K. Also the situation is not very different in rotating 
and non rotating models, the central temperatures and densities 
being weakly affected by rotation. 

This study emphasizes the importance of some reactions: 
proton-captures on 20 Ne and 23 Na respectively govern the 
amount of sodium produced and the linkage between the NeNa 
and MgAl chains. Concerning the MgAl chain, the key reaction 
is the proton-capture on 24 Mg. A raise of its rate with respect to 
the published values is required to reproduce the extreme vari- 
ations of magnesium and aluminum in the central region of a 
60 M star. This is due to the fact that the required temperature 
(of about 72-78 xl0 6 K, see Prantzos & Charbonnel in prepara- 
tion) is reached in the stellar core only at the very end of the 
main sequence in this object (see Fig.[7]l. 

4. Mixing and ejection of matter in a rotating 60 M Q 
star 

In the previous section we discussed what happens from the nu- 
cleosynthetic point of view in the core of a (rotating or non- 
rotating) 60 M star. We will now describe the critical effects 
of rotation on stellar mass loss and internal chemical structure. 
Finally we will follow the evolution of the surface abundances 
and compare our predictions with the data in GC stars. 

4. 1 . Mass loss 

Fig. [4] shows the evolution with time of the ratio between the 
surface and the break-up velocity in our 60 M Q model with 
Vi„i = 800 km s _1 . This ratio lies around 0.95 at the begin- 
ning of the main sequence. Rapidly, it falls down to 0.91 due 
to the establishment of differential rotation in the star: during 
this initial phase angular momentum is carried from the ex- 
tern al regions into the intern al layers and the core accelerates 
(see iMevnet & MaederlfeOOOl) . Later on Q/O cl it increases up to 




0.6 0.4 0.2 

x c 

Fig. 4. (top) Ratio of the surface rotation velocity to the break- 
up value for model 60rC. (bottom) Evolution of the total stellar 
mass for models 60rC and 60sC (full and dotted lines respec- 
tively) along the main sequence (the abscissa is the central mass 
fraction of hydrogen) 



reach unity (see § 2). However we force the model to stay be- 
low this level in order to avoid numerical difficulties (see § 2.2). 
Fig.|4]shows that all the episodes of strong mass loss correspond 
to a decrease of the surface velocity below the critical value. 
Then surface rotation velocity raises again. This process proba- 
bly leads to the formation of a circumstellar disk. We make the 
hypothesis that this material is eventually lost by the star. 

During the whole main sequence mass loss is dominated by 
this strong rotation-induced mechanical wind and is increased by 
more than a factor of 24 with respect to the standard case (e.g., 
model 60rC loses 20 M whereas the non rotating model expels 
less than 1 M ). 

As we just described, fast rotation does change the stellar 
mass loss quantitatively. But it has also a crucial qualitative im- 
pact on the properties and the topology of the ejecta. Indeed 
when the star rotates close to or at break-up the centrifugal force 
balances gravity. As a result, the mass loss is at least partly 
mechanically-driven, and the equatorial matter is released into a 
Keplerian disk. Very likely an equatorial disk forms as observed 
around Be stars (e.g JPorter & Riv inius 2003). This ejected ma- 
terial will thus be very easily retained within the GC potential 
well. This is very different from the non-rotating situation where 
the radiative winds escape at high velocity. This is one of the key 
points of our scenario which will be discussed in detail in § 7. 

As explained in § 2.2 the model evolves away from the Q.- 
Limit when it leaves the MS. At that moment however the most 
massive stars (M > 60 M ) encounter the QF-Limit. In the 
present work, we suppose that the matter lost at the Qr-Limit 
is released in an equatorial disk, as in the case of the material 
lost at the Q-Limit. When the star moves away from the QT- 
Limit due to heavy mass loss, the radiatively-driven fast winds 
take over. 
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Fig. 5. Abundance profiles at the end of main sequence (X c = 0.01) for the 60 M models in the standard case (model 60sC, top) 
and for V ini = 800 km s _1 (model 60rC, bottom) 



4.2. Abundance profiles 

Since mass loss is very efficient in the rotating case, the star peels 
off revealing the H-processed layers at the surface. Rotational 
mixing strengthens the modifications of the surface abundances 
as the products of central burning are transported outwards by 
various instabilities induced by rotation|f| 

In Fig.|5]we show the abundance profiles of various chemical 
elements as a function of the Lagrangian mass for the 60 M Q star 
(computed using set C for the nuclear network) at the end of the 
main sequence in the standard and rotating models (Vj n ; = and 
800 km s -1 respectively). In both cases the size of the convec- 
tive core is very similar (respectively 23.7 and 23.5 M in the 
standard and rotating models at the end of the main sequence) 
as well as the final central abundances (the central temperature 
differs by less than 1% between the two models). 

Outside the convective core however the abundance profiles 
show some notable differences. The most striking effect is due 
to mass loss which is strongly enhanced in the rotating model. 
The induced strong "peeling" exposes the nuclear regions. As a 
consequence the surface abundances and wind at the end of the 
main sequence hold the signatures of H-processing. 

4.3. Composition of the ejecta 

During the evolution of our 60 M rotating models we can dis- 
tinguish the following phases: 



5 The most important ones are the mer idional currents and the shear 
turbulence (see the review by |Talorl2004l) . The inclusion of these trans- 
port mecha nisms has improved the massive star models in many re- 
spects (see Maeder & Meynet 2 0001 for references); in particular, ro- 
tating models can re produce the chemical enrichment observed at th e 
surface of OBA stars dHeger & L anger 2000l: lMevnet & Maedej|200fjh . 



1 . The very beginning of the main sequence during which the 
surface velocity approaches the critical limit (between X c ~ 
0.76 to 0.65 in Fig. H. 

2. The rest of the main sequence when the surface velocity is at 
(or very near) the critical limit (for X c below 0.65 in Fig. |4j». 

3. The beginning of the core He-burning phase when the stellar 
luminosity is close to the Eddington value and the surface 
velocity stays close to the Or-Limit. 

4. The end of the He-burning phase. From this moment on the 
star is away from the Qr-Limit and the mass loss is mainly 
radiatively-driven. 

During the phases 2 and 3 described above, mass loss is mainly 
mechanically-driven and the stellar winds are supposed to be 
slow; they are enriched in H-burning products only. When the 
phase 4 starts the stellar surface and the ejecta are enriched in 
both H-burning and He-burning products, the winds are fast. 

Fig. [6] presents the evolution of the surface abundances in 
the rotating 60 M models as a function of the remaining stellar 
mass from the zero age main sequence up to the end of cen- 
tral He-burning; the different lines correspond to models com- 
puted with the various sets of nuclear reactions described in § 2. 1 
(60rA, B, C and D). We also indicate the range of abundances 
exhibited by the low-mass stars in NGC 6752 (shaded boxes). 

The surface abundance variations do mimic the central ones 
(see Fig. [3] and the discussion in § 3) with some delay due to 
non-instantaneous rotational mixing. In all the cases the surface 
abundances of carbon and oxygen at the end of the main se- 
quence fall respectively by 0.9 and 1 dex with respect to their 
original values, whereas that of nitrogen raises by 1.5 dex. The 
sodium surface abundance increases by 0.8 or 1.6 dex depending 
on the assumed nuclear reaction rates. Aluminum increases in all 
the cases but in a more efficient way when one enforces 24 Mg- 
burning by increasing significantly the corresponding reaction 
rate (set D). In summary and as discussed previously, the O-Na 
anticorrelation naturally shows up in all the cases. However one 
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Fig. 6. Evolution of the surface abundances as a function of re- 
maining mass for the 60 M Q rotating models computed with 
various sets of nuclear reactions : set A (full line), set B (long 
dashed lines), set C (short dashed lines) and set D (dotted lines). 
The short-dashed vertical line indicates the end of the main se- 
quence. The Or-Limit phase is located between the short- and 
long-dashed lines. The dotted line marks the moment when He- 
burning products start to appear at the surface. The shaded boxes 
on the left indicate the amplitudes of the abundance variations 
observed at the surface of low-mass stars in NGC 6752. 

has to call for an increase of the 24 Mg(p,y) reaction rate in order 
to correctly build up the Mg-Al anticorrelation. 

The H-burning signatures are reinforced at the stellar surface 
at the beginning of the central He-burning phase mainly as a 
result of previous mass loss episodes. 

However later on the products of He -burning start to show up 
at the surface. In particular the oxygen abundance raises steeply 
from 10~ 5 to nearly 0.6 in mass fraction when the total stellar 
mass is lower than ~ 28 M . C and N also display strong en- 
hancement at this phase. After the abundance of nitrogen de- 
creases being transformed into 22 Ne. 

In summary the ejecta of fast rotating massive stars display 
the chemical patterns observed in GC stars both on the main se- 
quence and during the fiF-Limit phase. Importantly, these are 
exactly the phases where matter is ejected by gently blowing 
winds which can be easily retained within the GC potential well. 
At the opposite, the He-burning products are released in the in- 
terstellar matter through fast radiatively-driven winds. We will 
come back to that crucial point of the WFRMS scenario in the 
discussion. 

4.4. Effect of the initial rotation velocity 

In order to investigate the impact of the initial rotation veloc- 
ity on the overall predictions we have computed a 60 M model 
with an initial velocity of 600 km s _I (model 60rE, with set C) 
instead of 800 km s _1 (model 60rC). The mixing of chemicals is 
less efficient in the slowest model as the shear turbulence, which 
dominates the transport processes, decreases with lower angu- 



lar velocity. However model 60rE reaches the break-up velocity 
later during the main sequence (when the central hydrogen mass 
fraction is 0.4 instead of 0.65 in model 60rC). Both effects tend 
to compensate each other: the mixing has a smaller efficiency but 
has more time to act before the star reaches the break-up. This 
leads to abundance variations which are even more pronounced. 
As an example the nitrogen abundance in winds when break-up 
is reached is 55% higher in model 60rE than in model 60rC (see 
Table©. 

When the stars are at break-up, the mass loss rate is higher 
for a lower initial rotation velocity: compared to the standard 
(non-rotating) case it is 29 and 24 times higher in models 60rE 
and 60rC respectively. The total stellar mass at the end of the 
main sequence is higher (by ~ 3 M ) in model 60rC. These dif- 
ferences play only a minor role for the composition of the yields 
which differ by only a few percents between the two models. 

Let us note finally that with an initial velocity around 
300 km s _1 the 60 M Q star would fail to reach the break-up 
and the conditions for the WFRMS scenario would not be ful- 
filled. If we infer models with various initial mass and metallicity 
(Ekstrom et al, in preparation), the 60 M would reach the criti- 
cal velocity with initial velocity as low as 400 km s . However 
fast rotation is possibly a characteristic of massive stars in a very 
dense environment such as forming GC where multiplicity and 
stellar collisions can play an important role. This point is dis- 
cussed in more detail in § 7. 

4.5. Summary 

To sum-up, the WFRMS appears to be very promising. 
Rotational mixing efficiently transports elements from the con- 
vective core to the surface. A high initial rotation velocity al- 
lows the star to reach break-up early on the main sequence and 
to eject important quantities of material loaded with H-burning 
products. This material is probably ejected through a slow wind 
and has great chance to remain in the GC potential well (see § 7). 

The abundance patterns at the stellar surface and in the ejecta 
follow those created in the core with some delay. The amplitude 
of the predicted O-Na anticorrelation reproduces well the ob- 
servational feature. Although the 60 M star produces too few 
Al when one sticks to the nuclear experimental value for the 
24 Mg(p,y) reaction rate, an increase of this rate allows one to 
explain the observed Mg-Al anticorrelation. 

5. Dependence on the initial stellar mass 

Let us now discuss how the theoretical predictions depend on 
the initial stellar mass. Standard and rotating models with initial 
masses ranging between 20 and 120 M are presented; in all the 
cases we used the set C for the nuclear reaction rates. 

5. 1 . Nucleosynthesis and mixing 

The initial mass has a direct effect on nucleosynthesis through 
the changes in the central temperature (see Fig.[7]l. As a result the 
NeNa and MgAl chains are more active in the warmer convective 
core of more massive stars. 

The evolution along the main sequence of the central abun- 
dances of some interesting nuclei is shown in Fig. [8] for the 20 
and 120 M . It can be compared to Fig. [3] The following differ- 
ences can be noted: 

- The CNO equilibrium value of O is slightly lower in the 
higher temperature regime of the 120 M stellar model. 
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Fig. 7. Evolution of the central temperature for the 20, 40, 60, 
120 and 200 M rotating stars as a function of the central hydro- 
gen abundance during the main sequence. 
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Fig. 8. Evolution of the central abundances of 16 0, 23 Na, 24 Mg 
and 27 Al as a function of the central hydrogen abundance on the 
main sequence for the 20 and 120 M rotating models. 



- 23 Na is first produced to the same extend in both models, but 
later on it decreases faster in the more massive star. 

- In the 120 M model 24 Mg decreases by 1.2 dex at the very 
end of the H-burning phase, while the 20 M never reaches 
sufficiently high temperatures for this element to burn. 



Rotating models 




50 100 50 100 



M ini (M Q ) M ini (M ) 

Fig. 9. Composition of the matter ejected in the winds as func- 
tion of initial stellar mass for the rotating and standard models 
(upper and lower panels respectively). The ejecta are presented 
in solar mass (left) or scaled with respect to initial mass (right). 
The winds enriched in H- and He-burning products are quoted 
respectively as "H wind" and "He wind". "Origin" refers to mat- 
ter which is not modified by nuclear burning. "Final mass" indi- 
cates the mass of the stars at the end of central He-burning. 

- 27 Al is produced earlier in the 120 M model. Nevertheless 
the plateau stays at the same level in both models. Only at 
the very end of the main sequence, does the Al abundance 
rapidly increase in the more massive star as a result of 24 Mg 
burning. 

Thus the main relevant difference between both stars con- 
cerns the MgAl chain. Only in the 120 M model and at the 
very end of central H-burning are 24 Mg and 27 Al respectively 
destroyed and produced. Even the 200 M star does not reach 
central temperature high enough to convert 24 Mg before the very 
end of main sequence. 

As discussed in lMevnet & Maederl fcOOd) . rotational mixing 
is more efficient when the initial stellar mass is higher, favoring 
the transport of the nuclear products outwards. At a given evo- 
lutionary stage, more massive stars also present stronger winds, 
both radiatively- and mechanically-induced. 

5.2. Winds 

Fig. [^illustrates the differences between the standard and rotat- 
ing models for the wind composition and the remaining stellar 
mass at the end of central He-burning. As already discussed, the 
rotating models lose much more mass than the standard ones 
due to three reasons: first, they do undergo mechanical winds 
when they reach the break-up velocity; second, their radiatively- 
driven winds are enhanced by the correcting factor due to rota- 
tion; third, they enter the WR phase earlier in their evolution. 

The rotating models with an initial mass above 40 M lose 
about half of their mass through winds loaded in H-burning 
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Fig. 10. (Top) Crosses are observed [O/Na] r atio in individ- 
ual NG C 6 752 stars computed from the data of iGrundahl et alj 
(120021) and lYong et al.l d2003l) . (Other panels) Theoretical his- 
tograms built from the wind composition of rotating massive 
stars in which 30% of pristine matter is added as required by 
the Li behavior (see the text). The white and shaded areas indi- 
cate respectively the slow (i.e., when the star is at breakup) and 
fast winds 

products. Moreover this matter is released very smoothly in the 
interstellar medium, with velocities low enough to remain in 
the GC. This is not the case of the winds of non-rotating mas- 
sive stars. Indeed in that case, the wind velocities on the main 
sequence are of th e order of 1000 to 2000 km s _1 (see e.g. 
lLamers et al. 1 119951) and this material is probably thrown off the 
gravitational potential well of the cluster. 

6. Comparison with observed abundance variations 
in GCs 

Here we explore how the ejecta of massive rotating stars can 
be used to form low-mass stars displaying abundance variations 
in light elements. It is out of the scope of the present paper to 
present a detailed study of the interaction between the disk en- 
closing fast rotating massive stars and the surrounding proto- 
cluster gas. We will thus only make some simple assumptions to 
verify if massive stars can be the progenitors of the long-lived 
stars we observe today. 

As the disks may be mixed with some pristine gas, we need 
first to evaluate the amount of this dilu tion. For that we use 
abundance variations in Li measured by iPasquini et al.l d2005l) 
in NGC 6752: Li is found to be as low as A(Li) = 1.93 (Alonso 
temperature scale) in the most polluted stars. In the present mod- 
els we do not follow explicitly the nucleosynthesis of lithium but 
we can reasonably assume that this fragile element is completely 
destroyed due to high temperature is the stellar interior and to the 
strong mixing induced by rotation. If we assume that the initial 
lithium abundance of the intracluster gas corresponds to the cos- 
mological value of 2.61 dCoc et al. 2004) and that no lithium is 



15 

10 

5 


30 

20 

10 

-§ o 




10 



NGC 675< 







1 1 1 1 


60rD 








II 



120rC 



I I I 





60rC 








1 1 1 1 

















40rC 


MM 


MM 


1 1 1 1 




Fig. 1 1. Same as Fig.lTOlfor 
120051) 



-1 -0.5 0.5 

[C/N] 

[C/N] ratio (data from lCarretta et al.l 





6 




4 




2 









20 




10 







S 


10 
8 




6 


■o 

C 


4 


i 


2 


S 







8 
6 




4 




2 









3 




2 




1 








NGC | 6752 | 



a 



60rD 



+ 



-H- 



120rC 
I | M I 

60rC 




-0.5 
[Mg/Al] 



Fig. 12. Same as Fig. [TOl for fM g/All ratio (data from 
IGrundahl et ai1l2002l: lYong et al.ll2Q03h 



present in the wind of massive stars, we need to add about 30% 
of pristine gas to those winds in order to find the extreme value 
of lithium (see also Prantzos & Charbonnel in preparation). 

As a first estimate we use this dilution factor between the 
stellar winds and the original matter and check whether it allows 
to reproduce all the observed abundance variations. In real proto- 
globular cluster there will surely be a dispersion in the amount 
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of dilution which in turn creates a variation in the abundances of 
low-mass stars. 

Fig. [10] [TTJ and [12] display the composition of the mixed 
matter formed by adding 30% of pristine gas to massive star 
ejecta as well the observed abundance dist r ibutions in NGC 675 2 
(iGrundahl et al.l l200l lYong et al.1 l200l ICarretta et al.l 120051) . 
The theoretical histograms indicate the mass ejected by winds 
mixed with pristine gas resulting on a given chemical composi- 
tion. The white areas correspond to the mass lost at break-up, 
while the hatched regions indicate that the corresponding winds 
are released when stellar rotation is well below its critical value 
(this occurs at the very beginning of the main sequence and dur- 
ing the He-burning phase). This latter component has thus a high 
velocity and probably escapes the GC potential well. Only the 
H-rich loaded winds are shown here. The He-rich winds are re- 
leased at high velocity and are supposed to escape the potential 
well of the GC. 

Let us first concentrate on the [O/Na] distribution. The abun- 
dance pattern in the mixed matter of the 60 and 120 M mod- 
els covers the entire observational range. The case of the less 
massive stars is different. At first glance the predictions for the 
40 M model also account well for the observed dispersion. This 
is not the case when one considers only the slow winds ejected 
at break-up. The 20 M Q star follows the same trend but with a 
smaller extent. Let us note that the highest [O/Na] ratio for both 
the 20 and 40 M stars is shifted to the left compared to more 
massive models. This is due to the longer time required to reach 
critical velocity at the beginning of the main sequence. When the 
break-up is reached rotational mixing had more time to transport 
the elements through the radiative envelope. 

Turning our attention to the [C/N] distribution (Fig. [TTJ we 
find that the ejecta of rotating massive stars cover very nicely the 
observational range. Again the 20 M and 40 M , being more 
mixed before reaching break-up, do not expel matter with origi- 
nal composition (see Table[5]). 

Regarding the [Mg/Al] ratio we see that the winds of all the 
models computed with the set C (experimental limits) present 
lower abundance variations than required by the observed distri- 
bution. In addition let us recall that the magnesium isotopic ra- 
tios obtained in these models are at odd with the observed ones: 
even in the case of the 120 M aluminum is built mainly from 
25 Mg and 26 Mg. Only model 60rD covers the whole observa- 
tional range both in terms of Mg and Al abundances and of Mg 
isotopic ratios. The same would be true for other masses com- 
puted with set D. This difficulty remains even in the extreme 
case, where present day halo stars would be naked cores of mas- 
sive stars having undergone strong evaporation at the end of the 
main sequence. Only at the end of main sequence the 60 and 
120 M models manage to destroy 24 Mg by more than 0.3 dex 
and to enhance 27 Al. However at that time the high central tem- 
perature leads to a strong destruction of 23 Na which falls near or 
under its initial value. 



7. Discussion 

Let us now discuss some details of the WFRMS scenario. 

7.1. Kinematics and topology of the ejected material 

Figure [13] shows the distribution of the globular clusters as a 
function of the escape velocity estimated at the center and at 
half-mass radius of t he present-day cluste r (continuous and dot- 
ted line respectively. iGnedin et al.l [20021) . One sees that more 
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Fig. 13. Histogram of the number of globular clusters having 
a given value of th e escape velocity. The data are taken from 
IGnedin et al.l d2002l) . The continuous line shows the results for 
the escape velocity at the center of the cluster, the dotted line at 
the cluster half-mass radius. 



than 60% of the 141 globular clusters used to build the histogram 
(i.e. the large majority of the known galactic globular clusters) 
have an escape velocity at the center between 20 and 1 00 km s~ 1 . 
The range is shifted to lower values if escape velocities estimated 
at half-mass radius are used. Thus, taken at face, these values 
imply that all the wind ejecta of O-type stars and of WR stars 
(wind velocities of a few thousandths of km s~'), of Luminous 
Blue Variable (wind velocities of a few hundred km s~') as well 
as the supernova ejecta (velocities of the ejecta of the order of 
10000 km s~') will be lost by the globular cluster. This state- 
ment must however be a little tempered by a couple of facts: first 
the escape velocity of the present day globular clusters may be 
different from the one at the time of their formation. Would the 
present globular clusters be the remnants of much more mas- 
sive systems having lost a great part of their initial mass, then 
the present day escape velocity might be a poor estimate of the 
escape velocity at the time of their formation^ Second, the inter- 
stellar medium at the time of the globular cluster formation was 
much more dense and thus the matter ejected was slowed down 
by the bow shock between the ejecta and the a mbient medium as 
shown for instance in the numerical models bv lFrever et al.l2003l 
(see also lGarcia-Segura et al.ll996h . It is however difficult, with- 
out a more precise view of the initial conditions of the globular 
cluster formation to come to firm conclusion on these pointsQln 



6 More precisely we need here the escape velocity at the time of the 
chemical inhomogeneities formation, which in our present model would 
be different from the time of formation of the first stellar generation. 

7 Of course one could argue that we can observe young globular 
clusters and thus obtain some constraints on their formation. However, 
many (if not all) young super star clusters observed today and believed 
to be young globular clusters, result from galaxy merging or interac- 
tions. It is not clear at all that the old galactic globular clusters were 
formed in this way. 
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Table 5. Sum-up of the main features of our models: initial mass, initial surface velocity, and for each evolutionary phases the 
remaining mass and the mean mass fraction in the wind for several isotopes integrated from the zero age main sequence. The 20, 
40 and 120 M models are computed with the reaction rates from set C. For rotating models the last division indicates abundance 
when He-burning products appear at the surface, or at the end of He-burning phase depending on the first case reached. 
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Standard Models 
























End of central H-burning 












20 





19.9 


7.5e-l 


2.4e-l 


3.5e-5 


1.0e-5 


3.0e-4 


1.5e-8 


3.3e-7 


1.7e-5 


2.1e-6 


2.1e-6 


9.0e-7 


40 





39.7 


7.5e-l 


2.4e-l 


3.5e-5 


1.0e-5 


3.0e-4 


1.5e-8 


3.3e-7 


1.7e-5 


2.1e-6 


2.1e-6 


9.0e-7 


60 





59.2 


7.5e-l 


2.4e-l 


3.5e-5 


1.0e-5 


3.0e-4 


1.5e-8 


3.3e-7 


1.7e-5 


2.1e-6 


2.1e-6 


9.0e-7 


120 





117.0 


7.5e-l 


2.4e-l 


3.5e-5 


1.0e-5 


3.0e-4 


1.5e-8 


3.3e-7 


1.7e-5 


2.1e-6 


2.1e-6 


9.0e-7 












End of central He-burning 












20 





19.9 


7.5e-l 


2.4e-l 


3.5e-5 


1.0e-5 


3.0e-4 


1.5e-8 


3.3e-7 


1.7e-5 


2.1e-6 


2.1e-6 


9.0e-7 


40 





39.1 


7.5e-l 


2.5e-l 


3.5e-5 


1.0e-5 


3.0e-4 


1.5e-8 


3.3e-7 


1.7e-5 


2.1e-6 


2.1e-6 


9.0e-7 


60 





45.9 


6.9e-l 


3.1e-l 


2.6e-5 


7.9e-5 


2.3e-4 


1.2e-8 


1.2e-6 


1.7e-5 


1.6e-6 


1.7e-6 


1.8e-6 


120 





89.3 


5.9e-l 


4.1e-l 


1.5e-5 


1.7e-4 


1.4e-4 


6.3e-9 


2.0e-6 


1.7e-5 


9.9e-7 


1.2e-6 


3.0e-6 














Rotating Models 


























Break-i 


up reached 














20 


600 


20.0 


7.5e-l 


2.4e-l 


l.le-5 


4.9e-5 


2.9e-4 


1.5e-8 


7.5e-7 


1.7e-5 


2.1e-6 


2.2e-6 


9.3e-7 


40 


800 


40.0 


7.5e-l 


2.4e-l 


3.5e-5 


1.3e-5 


3.0e-4 


1.6e-8 


3.7e-7 


1.7e-5 


2.1e-6 


2.2e-6 


9.2e-7 


60rE 


600 


59.7 


7.5e-l 


2.5e-l 


3.3e-5 


1.7e-5 


3.0e-4 


1.5e-8 


4 2e-7 


1.7e-5 


2.1e-6 


2.1e-6 


9.2e-7 


60rA 


800 


59.9 


7.5e-l 


2.4e-l 


3.5e-5 


l.le-5 


3.0e-4 


1.6e-8 


j.jc- / 


1.7e-5 


2.1e-6 


2.1e-6 


9.2e-7 


60rB 


800 


59.9 


7.5e-l 


2.5e-l 


3.5e-5 


l.le-5 


3.0e-4 


1.6e-8 


3.5e-7 


1.7e-5 


2.1e-6 


2.2e-6 


9.2e-7 


60rC 


800 


59.9 


7.5e-l 


2.4e-l 


3.5e-5 


l.le-5 


3.0e-4 


1.6e-8 


3.4e-7 


1.7e-5 


2.1e-6 


2.2e-6 


9.2e-7 


60rD 


800 


59.9 


7.5e-l 


2.5e-l 


3.5e-5 


l.le-5 


3.0e-4 


1.6e-8 


3.4e-7 


1.7e-5 


2.1e-6 


2.2e-6 


9.2e-7 


120 


800 


119.5 


7.5e-l 


2.4e-l 


3.5e-5 


l.le-5 


3.0e-4 


1.6e-8 


3.4e-7 


1.7e-5 


2.1e-6 


2.1e-6 


9.2e-7 












End of central H-burning 












20 


600 


18.4 


7.3e-l 


2.7e-l 


1.6e-6 


8.8e-5 


2.5e-4 


1.3e-8 


1.5e-6 


1.7e-5 


1.8e-6 


2.0e-6 


1.3e-6 


40 


800 


28.8 


7.2e-l 


2.8e-l 


1.3e-5 


9.3e-5 


2.3e-4 


1.2e-8 


1.4e-6 


1.7e-5 


1.7e-6 


1.9e-6 


1.6e-6 


60rE 


600 


42.7 


6.7e-l 


3.3e-l 


1.9e-5 


l.le-4 


2.1e-4 


1.0e-8 


1.5e-6 


1.7e-5 


1.5e-6 


1.6e-6 


2.0e-6 


60rA 


800 


39.9 


6.8e-l 


3.2e-l 


1.9e-5 


1.2e-4 


2.0e-4 


1.0e-8 


7.1e-6 


1.9e-5 


1.4e-6 


2.2e-6 


1.5e-6 


60rB 


800 


39.9 


6.8e-l 


3.2e-l 


1.9e-5 


1.2e-4 


2.0e-4 


1.0e-8 


6.4e-6 


1.7e-5 


1.4e-6 


2.0e-6 


1.7e-6 


60rC 


800 


39.7 


6.8e-l 


3.2e-l 


1.9e-5 


1.2e-4 


2.0e-4 


1.0e-8 


1.7e-6 


1.7e-5 


1.4e-6 


1.6e-6 


2.1e-6 


60rD 


800 


39.9 


6.8e-l 


3.2e-l 


1.9e-5 


1.2e-4 


2.0e-4 


1.0e-8 


1.7e-6 


1.6e-5 


1.5e-6 


1.8e-6 


2.7e-6 


120 


800 


76.6 


6.1e-l 


3.9e-l 


1.7e-5 


1.8e-4 


1.3e-4 


6.3e-9 


2.2e-6 


1.7e-5 


9.2e-7 


l.le-6 


3.1e-6 








Appearence of He- 


■burning 


products 


at the surface/End of 


central He-burning 








20 


600 


16.4 


7.0e-l 


3.0e-l 


1.6e-6 


l.le-4 


2.3e-4 


l.le-8 


1.8e-6 


1.7e-5 


1.7e-6 


1.9e-6 


1.6e-6 


40 


800 


18.5 


5.7e-l 


4.3e-l 


9.3e-6 


1.6e-4 


1.6e-4 


8.0e-9 


2.1e-6 


1.7e-5 


1.2e-6 


1.3e-6 


2.7e-6 


60rE 


600 


33.2 


5.2e-l 


4.8e-l 


1.4e-5 


1.7e-4 


1.4e-4 


7.2e-9 


2.0e-6 


1.7e-5 


l.le-6 


1.2e-6 


3.0e-6 


60rA 


800 


30.8 


5.4e-l 


4.6e-l 


1.5e-5 


1.7e-4 


1.5e-4 


7.2e-9 


8.6e-6 


2.5e-5 


1.3e-6 


2.0e-6 


2.5e-6 


60rB 


800 


30.7 


5.5e-l 


4.5e-l 


1.5e-5 


1.7e-4 


1.5e-4 


7.2e-9 


1.0e-5 


1.7e-5 


l.le-6 


1.5e-6 


2.6e-6 


60rC 


800 


30.7 


5.5e-l 


4.5e-l 


1.5e-5 


1.7e-4 


1.5e-4 


7.3e-9 


2.0e-6 


1.7e-5 


l.le-6 


1.2e-6 


3.0e-6 


60rD 


800 


31.0 


5.5e-l 


4.5e-l 


1.5e-5 


1.7e-4 


1.5e-4 


7.3e-9 


2.0e-6 


1.4e-5 


l.le-6 


1.6e-6 


5.5e-6 


120 


800 


68.7 


5.0e-l 


5.0e-l 


1.5e-5 


2.0e-4 


l.le-4 


5.2e-9 


2.2e-6 


1.7e-5 


8.4e-7 


9.6e-7 


3.5e-6 



the following, in absence of a better knowledge of the globular 
cluster formation conditions, we shall consider that the present 
escape velocities are representative of the escape velocities at the 
time of the chemical inhomogeneities formation. 

Thus we have to look if our massive star models are able to 
eject matter at sufficiently low velocities for allowing the ejecta 
to remain in the potential well of the globular cluster. Rotation 
might be of great help in this respect. Indeed, in situations where 
the surface equatorial velocity is such that the centrifugal ac- 
celeration exactly counterbalances the gravity, matter can be 
launched into a Keplerian orbit and an equatorial disk easily 
forms as schematically shown in Fig. [14] The material mainly 
"mechanically" ejected into the disc will be possibly available 
for further star formation, unless it falls back onto the star and/or 
is ejected outside the globular cluster by further violent wind 
episodes or by the supernova explosion. Let us briefly discuss 
these different possibilities: 

- Fall back of the disk mat erial? Be stars are main- sequence 
stars surrounded by a disc dPorter & Riviniusl2003l) . The ori- 



gin of the disk is probably stellar rotation near the critica l 
limit as convincingly discussed by iTownsend et al.l d2004l) . 
Thus these objects are exactly the observable counterparts 
of the stellar mod els we are inter ested in here. Observations 
(see the review byli vimu 120051) indicate that Be star disks 
are eroded and finally dissipate, adding their material to the 
ambient interstellar medium. Thus these observations would 
favor the lost of the disk material around Be stars and not its 
fall back onto the star. 
- Ejection by fast stellar winds? When the star does no 
longer rotate with the critical velocity, matter is no longer 
preferentially ejected in the equatorial plane, and different 
situations may occur. First it might be that the timescale for 
disk dissipation is very short and thus that most of the disk 
material has already disappeared when the faster radiatively 
driven winds set in. In that case there is the possibility that 
the fast wind pushes this freshly ejected material out of the 
globular cluster, which will be then lost. The most favorable 
case for keeping the material inside the globular cluster is 
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Fig. 14. Schematic view of the WFRMS scenario showing the 
possible geometry of the stellar ejecta at various evolutionary 
phases: a) during the main sequence, when the star rotates near 
or at the critical velocity, matter is preferentially ejected in the 
equatorial plane by the action of the centrifugal acceleration; b) 
after the main sequence, the surface velocity is no long criti- 
cal and the wind is triggered mainly by radiation. It is no long 
equatorial and becomes isotropic; c) the supernova explosion re- 
sulting from an initially fast spinning star, if it occurs, may favor 
ejection through jets aligned along the rotational axis. 

that the disk has a sufficiently long lifetime for still keep- 
ing most of its material when the fast winds set in. In that 
case, if the star is still in the blue part of the HR diagram 
(as is the case for our 60 and 120 M ) and if its rotation is 
still sufficiently close to the critical value, the wind will be 
preferentially ejected along the rotational axis. This would 
thus prevent the winds to sweep off the disk. If the star has 
a rotational velocity well below the critical one, the wind 
is isotropic and one may hope that the fraction of the wind 
which may interact with the disk will be too small to destroy 
it (see Fig. IT4T>. 

- Will the disk be thrown away by the supernova explo- 
sion? In case a black hole is formed that swallows all the 
mass of the pre-supernova, no explosion occurs and the 
disk will not be affected. If, on the other hand, a super- 
nova explosion occurs, the fast rotation of the core may fa- 
vor ejection along the ro tational axis as in the models of 
iMaeda & Nomotd ( l2003l) . In that case there is some chance 
that the disk will not be affected by the supernova explosion 
as illustrated schematically in Fig. [14] Of course the disk ma- 
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Fig. 15. Evolution of the ratio between the surface angular ve- 
locity and the classical critical velocity (top) and of the surface 
equatorial velocity (bottom) as a function of the remaining mass 
for the 60 M star. Vertical short-dashed line corresponds to the 
end of the main sequence and dotted line indicates the moment 
when He-burning products show up at the surface. The fir- 
Limit phase is delimited by the short- and long-dashed lines. 



terial can also be destroyed and ejected out of the GC pushed 
away by the winds and/or supernova shocks from nearby 
massive stars. Our scenario will then be possible only if in 
relatively short timescales (shorter than typically the aver- 
aged time between two supernova events), new stars form in 
the disk, or more probably in the interaction region between 
the disk material and the interstellar medium. In that case, the 
fast moving material will escape following the nearly empty 
channels left behind by the star formation process, in line 
with the scenario proposed by PC06. 

These ideas are highly speculative and will have to be checked 
by detailed hydrodynamical calculations in the future. For the 
purpose of the present discussion, we shall suppose that the 
equatorial disk formed when the star reaches the critical limit 
consists of material which will be made available for forming 
new stars in the globular cluster. 

Now, in the frame of the above hypothesis, we have to see 
when this occurs during the evolution of our stellar models and 
whether sufficient mass with the adequate chemical composition 
is ejected during these phases. The present models show that the 
critical limit can be reached at two epochs: 

1 . During the main sequence, provided that the initial velocity 
is high enough. From Fig.[T31one sees that our 60 M model 
with v m i = 800 km s _1 remains near the break-up limit dur- 
ing nearly the whole main sequence. At the turnoff, the star 
has lost more than 20 M Q . 
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2. After the main sequence we saw that stars massive enough 
do reach the QF-limit which maintains them at the break-up0 
due to the raise of the luminosity. Strong mass loss ensues: 
During this phase the 60 M star expels more than 5 M Q 
(see Fig.[T5l). We suppose here that this material will still be 
preferentially ejected along the equatorial plane and will join 
the disk of slowly outflowing material. 

7.2. Composition of the ejected material 

In the previous sections we have seen that the chemical compo- 
sition of the ejected material presents many similarities with the 
abundance patterns observed at the surface of globular cluster 
stars. If one focuses on the matter ejected during the main se- 
quence, one sees however that the theoretical variations can ac- 
count only partly for the observed range in oxygen and sodium. 
However slightly later in the evolution, strongly CNO processed 
material is ejected and the most extreme observational cases can 
be reproduced (see Table[5]l. The observed variations of Mg and 
Al can also be reproduced, provided the rate of the 24 Mg(p, y) 
reaction is enhanced by about a factor 1000 for temperatures 
around 50 x 10 6 K. However if the non-enhanced rate for this 
reaction is correct, then our models reproduce only part of the 
Mg-Al anticorrelation. 

One has to conclude that fast rotating massive stars are ex- 
tremely good candidates for providing the material from which 
the long-lived low-mass stars that we are presently observing 
formed. Of course the process leading to the incorporation of 
this material in new stars is far from being described here, but 
at least we have given some reasons to believe that the interest- 
ing ejecta can be retained in the globular cluster. It has also to 
be stressed that the stars forming from the wind ejecta of mas- 
sive stars mixed with some pristine interstellar gas wo uld also 
present higher helium abundances (ISalaris et al.l I2006D and of 
course high C/N ratios. The values of the 12 C/ 13 C would be quite 
low: all the massive stars have ratio below 5. The corresponding 
implications will be studied in a subsequent work. 

7.3. Why only in GCs? 

To end this section, we can wonder why such a process did only 
occur in GCs. The environment of field halo stars is of course 
very different from the one of the globular clusters. However, in 
the past, they probably share, at least for a short time, a cluster 
environment. These clusters were probably much less massive 
and/or much less concentrated than the progenitors of the present 
day globular clusters. This is required for these clusters to have 
evaporated, being completely disrupted either by the energy in- 
jections of the first supernovaefl or by tidal effects. In that case, 
star formation triggered by the wind/SN shocks from evolving 
massive stars might have been quenched, rendering impossible 
the birth of field stars with processed material. 

Another possibility, also related to the high stellar density 
in GCs, could be that such an environment favors higher rota- 
tional velocities than less dense environments. It is however not 
obvious that denser environments would favor high rotational 
velocities, the stellar encounters being able of both spinning up 
the stars or slowing them down. The only thing we can say at 



8 At this stage, the critical velocity is lower ed with respect to its clas- 
sical expression (see lMaeder & Mevnetll2000h . 

9 It is interesting to note that the gravitational energy of a globular 
cluster is of the order of 10 51 erg i.e. of the same order of magnitude as 
the kinetic energy emitted by a core collapse supernova explosion. 



the moment is that some observation s find higher ro t ationa l ve- 
locities in clusters than in the field: iHuang &~G ies ( 200fJ) de- 
termined the projected rotational velocities of 496 OB stars in 
clusters within the approximate age range 6-73 Myr. They found 
that there are fewer slow rotators amon g the cluster B-type stars 
relative to nearby B stars in the field. IStrom et all (|2005) also 
found that stars in h and x Persei tend to rotate faster than the 
fie ld stars counterpart s. A similar conclusion has been reached 
bv lDuftonetail ( l2006l) on the basis of the rotational velocities of 
stars in the two clusters NGC 3293 and NGC 4755. Of course it 
is difficult at the present time to draw firm conclusions about the 
origin of this difference. Let us just note that these observations 
may support the view according to which stars born in dense en- 
vironments may have faster rotational velocities than stars born 
in lose aggregates. 

An alternative to this view is to consider that all the field stars 
were born in the progenitors of the present day globular clus- 
ters. Let us recall that field halo sta rs contain at le ast hundred 
times the mass in globular clusters dWoltier|[l975l) . The above 
scenario would therefore imply a very efficient mechanism that 
would remove 99% of stars initially formed in the progenitor of 
the present day GC. Moreover, this mechanism should either be 
more rapid than the enrichment of the interstellar medium by 
massive stars, or should only remove stars with no chemical in- 
homogeneities. While it appears difficult to favor such a view, it 
would have the advantage of making the stars presenting inho- 
mogeneities a very small subset of the whole stellar population 
in a given cluster and thus alleviate the need of a very flat IMF 
(PCO6)0 

8. Conclusion 

The two main reasons why massive stars have been discarded in 
the past in the context of the GC self-enrichment scenario are 
related to the fact that these objects are expected to produce iron 
and to have fast winds. 

In the present paper, we propose the Wind of Fast Rotating 
Massive Star scenario in order to explain the chemical inhomo- 
geneities observed in GC stars. The key point of this scenario is 
the fast rotation of massive stars. Each one of these two char- 
acteristics is important: 

- Fast rotation is needed to remove material from the stellar 
surface and inject it with a low velocity in the interstellar 
medium. It triggers internal mixing which brings to the sur- 
face material processed in the core. This enables a star to 
eject material with chemical compositions similar to that ob- 
served at the surface of GC stars. 

- Massive stars have short lifetimes and can release H- 
synthesized material while low-mass stars are still forming 
in the nascent globular cluster. They can, through the wind 
and SN shocks or through the ionization front they produce, 
trigger star formation in their vicinity, being thus able to be 
at the same time the cause of new star formation and the 
provider of at least part of the material from which the stars 
form. Moreover, it appears that rotating massive stars can 
lose low-speed material only enriched in H-burning prod- 
ucts. Low- and intermediate-mass stars would not be able to 
do that, either because they are less efficiently mixed by ro- 
tation and/or they have more difficulty to reach the break-up 

10 Since massive stars are preferentially located in the central parts of 
a young cluster, they may have polluted mainly the central regions. The 
outer part of the GC progenitors would at least for a while be composed 
of non polluted stars which may have been stripped off by tidal effects. 
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limit and of course the Qr-Limit than the massive stars dur- 
ing the main sequence. In addition their central temperature 
is too low to efficiently activate NeNa and MgAl chains dur- 
ing the MS. Besides they require a more top-heavy IMF than 
massive stars as explained by PC06. 

Globular clusters may be suitable environments to form fast 
rotating massive stars in the proper range of metallicity. It is in- 
teresting to recall that the WFRMS scenario for explaining the 
GC abundance anomalies is only one consequence of fast rota- 
tion of massive stars. This kind of stars might be also interest- 
ing for understanding other features as the origin of the carbon- 
rich ultra metal-poor stars, the high N/O ratio observed in halo 
stars, and the high helium abundance of a part of stars in lo Cen 
dMaeder & Mevnetll2006l) . The anticorrelations observed in GC 
might thus be an additional observed consequences of the chem- 
ical enrichment expected from fast rotating massive stars. 
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